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Abstract 

We establish a rate of convergence of the two scale expansion (in the sense of homogenization theory) of 
the solution to a highly oscillatory elliptic partial differential equation with random coefhcients that are a 
perturbation of periodic coefficients. 



1 Introduction and presentation of the main result 

^^This article focuses on establishing a rate of convergence of the two scale expansion (in the sense of homoge- 
^inization theory) of the solution to a highly oscillatory partial differential equation with random coefficients. We 
CNbegin our exposition by briefly discussing the same question in a deterministic setting, before turning to the 
"^stochastic setting. 

Consider the highly oscillatory problem 

-div Aper (^-) Vuq = f in P, 

> , I ^^0 = ondV, 

• ^ - 

^where P is a regular bounded domain of M*^, / E L'^{'D), and Ap^r is a Q-periodic elliptic bounded matrix, with 
= (—1/2, 1/2)'^. For simplicity, we manipulate henceforth symmetric matrices, but the arguments carry over 
'to non-symmetric matrices up to slight modifications. It is well known (see e.g. the classical textbooks [7111^117]. 
and also [14J for a general, numerically oriented presentation) that Uq converges, weakly in (D) and strongly 
in L?'{V), to the solution Ug to 

-div [Al,,Vul\ = f in V, 

'1*0 = on dV, ^ ' 

where the homogenized matrix is given by 

= I {e^ + Vwl{y)fAper{y){ej + Vw^iy)) dy, (3) 
JQ 

where, for any p S W^, Wp is the unique (up to the addition of a constant) solution to the corrector problem 
associated to the periodic matrix viper : 



-div [Aper{p + Vwl)\ =0, 

)p is Q-periodic. 



(4) 
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The corrector function allows to compute the homogenized matrix, and it also allows to obtain a convergence 
result in the strong norm. Indeed, in dimension d > 1, under some regularity assumptions recalled below, 
we have 

'' d n 



Un 



dxi 



(5) 



m{T>) 



for a constant C independent of e (in dimension d 



Note that 



tin + e 



1, the difference is of order e rather than -v/e). 

At first sight, one could 



- ) — is a function of order 1 in the norm 



thus expect that the difference between Uq and Vq is of order e, rather than ^/e. This lower order (in dimension 
d > 1) is due to an inconsistency of the boundary conditions. Note indeed that, by definition, Uq = on dT>, 
which is not the case of Vq. Note also that the (lower than expected) rate in ([5]) is not specific to the choice of 
homogeneous Dirichlet boundary conditions in ([1]), and also holds for Neumann boundary conditions, as stated 
in da p. 29] (see also [22]). 

The order of approximation improves if we ignore the difference between Uq and Vq at the boundary of the 
domain (see [H Theorem 2.3]). Alternatively, one can build functions, the so-called boundary layers, that correct 

in the neighboorhood of dT), to eventually improve the accuracy of the approximation of Uq so obtained, 
in the complete domain V. We refer to [U EI] and to [131 Appendix B] (see also [21 Chap. 5] for the study 
of the same question in a time-dependent, parabolic setting). On another note, we refer to [26] for studies 
on the rate of convergence of Uq to Uq in the L°°{T>) norm (see also |14j and references therein, and [10] for 
extensions to some nonlinear cases), and to [2^122] for similar studies on the lowest eigenvalue Aq of the operator 



-div 



A. 



per 



The result ([5]) is interesting from the theoretical viewpoint. It is also helpful for proving numerical analysis re- 
sults. In particular, this result is a key ingredient to prove error bounds for the Multiscale Finite Element Method 
(MsFEM) . This numerical approach aims at approximating the solution Uq to the highly oscillatory problem ([1]) 
(for a small, but non vanishing small scale e), and does so by performing a variationnal approximation of ([1]) 
using pre-computed basis functions that are adapted to the problem. Consequently, the MsFEM approach yields 
an accurate approximation of Uq using only a limited number of degrees of freedom, in contrast to a standard 
Finite Element Method approach. In addition, the MsFEM approach is applicable in general situations, and 
is not limited to the case when the highly oscillatory coefficient of the equation reads A'^{x) = Ap^r ^ 

fixed periodic matrix j4per- See [13] and references therein. As described below, our motivation for this work 
stems from our work |19j . where we suggest a possible extension of the MsFEM approach to weakly stochastic 
settings. Again, a key ingredient for proving error bounds on the approach we propose there is to have a rate 
of convergence of the type 

Let us now turn to the stochastic case. As will be seen below, less precise results are known than in the 
deterministic, periodic case. The highly oscillatory problem reads 




w) Vu^(-,a;) 



/ 



in V, 
on dV, 



(6) 



where the matrix A^ is now a stationary symmetric matrix, uniformly elliptic and bounded (see ([5D below for 
a precise definition of stationarity, which is the common assumption in stochastic homogenization) . The role of 
the parameter r] will be made precise in ^ below. It can momentarily be ignored. Again, as in the periodic 
case, it is well known (see for instance [T7]) that converges, almost surely, weakly in H^CD) and strongly in 



2 



L?'{T>), to n*, solution to the homogenized equation 



-div [A^Vn^] = / 



in V, 
on &D, 



where the homogenized matrix is given by 



where, for any p G 
corrector problem 



\ Wp is the unique (up to the addition of a random constant) solution to the stochastic 

-div [Ar,{-,uj){p + Vw1{-,uj))] =0 in R'^, 
Vwp is stationary in the sense of ^ below, 



E 



Vw'^piy, •) dy 



0. 



As in the periodic case, the corrector function Wp allows to obtain a convergence result in the norm (see 
Theorem 3]): 



E 



u^(-,a;) 



£ 



du* 
dxi 



converges to as e — >• 0. 



(7) 



However, in contrast to the periodic case, the rate of convergence is generally not known, in dimensions higher 
than one. In the one-dimensional case, this question has been addressed in [9l |^. It is shown there that the 
rate can be arbitrary small, depending on the rate with which the correlations of the random coefficient in ([6]) 
vanish. The only assumptions of stationarity and ergodicity do not allow for a precise rate. See also [20j for 
the study of a similar question for a variant of stochastic homogenization, again in the one-dimensional case, 
and [5] for results in the multi-dimensional case, for a different equation. 

The aim of this article is to show that, in a weakly stochastic case (the precise sense of which is given 
below), a convergence rate for ([TD can be obtained (in the same spirit as ([5])). As in the deterministic case, 
this result is interesting from the theoretical viewpoint, and somewhat complements the one-dimensional results 
of [Hi O [20]. It is also useful from a numerical analysis viewpoint. In [19j, we propose an extension of the 
MsFEM approach to weakly stochastic settings, and we use there the homogenization result that we prove in 
this work (see Theorem [2] below) to obtain error bounds (see [191 Theorem 10]). 

Before presenting our result, let us briefly recall the basic setting of stochastic homogenization. Let (J7, J-", P) 
be a probability space. For a random variable X G L}{Q,,dF), we denote by E(X) = J^X{u})dF{u!) its expecta- 
tion value. We assume that the group acts on Q.. We denote by {Tk)k^id this action, and assume that 
it preserves the measure P, i.e. 

V/t G 1'^, G F, P(rfc^) = P(A). 

We assume that r is ergodic, that is, 

VA G T, (yk G Z^, TkA = a) ^ (P(A) = or 1). 



We define the following notion of stationarity: any F G L^^^ (M.'^, L^{Q)^ is said to be stationary if 

V/c G Z"', F{x + k,uj) = F{x,TkU}) almost everywhere, almost surely. (8) 
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Note that we have chosen to present the theory in a discrete stationary setting, which is more appropriate for 
our specific purpose, which is to consider a setting close to periodic homogenization. Random homogenization 
is more often presented in the continuous stationary setting. This is only a matter of small modifications. We 
refer to the bibliography for the latter. 

We now precisely describe the weakly stochastic setting we consider. We assume that the matrix Af^ in ^ 
reads 

Ar,{x,uj) = Aper{x) + r]Ai{x,uj), (9) 

where G M is small deterministic parameter, Aper is a symmetric uniformy elliptic bounded Q-periodic matrix, 
and Ai is a symmetric matrix, stationary in the sense of ([8|), and bounded: |74i(x,co')| < C almost everywhere 
in R*^, almost surely. We also assume that Aj^ is uniformly elliptic and bounded, in the sense that, for all 7/ G M, 

Ar,{-,uj) € {L^{R^)f''^ a.s. 

and there exists c > such that 

G M'', S,'^Ar,{x,uj)C > c C'^C a.s., a.e. on R'^. 
We furthermore assume that Ai is of the form 

Ai{x,Uj) = ^ lQ+kix)Xkiuj) Bper{x), (10) 

where {Xf^{uj))j^^^d is a sequence of i.i.d. scalar random variables such that 

3C, V/c e Z'^, \Xkiuj)\ < C almost surely, 

and Bper G (L~(M'^)) is a Q-periodic matrix. Finally, we assume that 

Aper is Holder continuous, (11) 
Bper is Holder continuous. (12) 

As pointed out above, the symmetry assumption is not essential, and our arguments below carry over to non- 
symmetric matrices up to slight modifications. Likewise, the assumption ()10p can be relaxed. What is important 
in ([To]) is that Ai is a sum of direct products of a function depending on x with a random variable, depending 
only on to. 

In contrast, it is difficult to weaken assumptions (jlip and ()12p . which are used to obtain some regularity 
on the correctors Wp and ipp, solutions to (jl]) and (fTUj) below, respectively. We indeed recall that, under 
assumption (jlip . we have G VF^'°°(R'^) for any peR'^ (see e.g. ^ Theorem 8.22 and Corollary 8.36]), and 
similarly for xpp, under assumption p^ . In the sequel, we will use the fact that Wp and ipp belong to W'^'°°{R'^), 
which is a standard assumption when proving convergence rates of two-scale expansions (see e.g. pj, p. 28]). 

We also note that, following [3], the assumption (|lip is useful to characterize the asymptotic behavior of the 
Green function associated to the operator L = —div [^perV-] on the domain D/e (with homogeneous Dirichlet 
boundary conditions). This Green function will be used in the sequel. 

Remark 1. There are several ways to formalize a notion of "weakly" stochastic setting, and Q is only one of 
them. We refer to fTB[ \^ for other examples. 

Our main result is the following. 
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Theorem 2. Assume that the dimension d is strictly higher than 1. Let be the solution to ([6]), and assume 
that Ar, satisfies Let A*^^, -w^ and Uq be defined by and il 

Ui € Hq{V) be defined by 



Let B £ 



pdxd 



yi <i,j<d, Bi 



/( 

JQ 



ei + VwlfBp^r{ej+Vwl 



and 



-div[A*^^Vu\\ = div[BVuQ\^ inV, 



= 



on dT>. 



and 
(13) 

(14) 



Introduce defined by 

d 



p=l 



+r?E(Xo)V'e, (-) dpu^o + 7/ E (Xkiuj) - E(Xo)) Xe, (- - k) d^u^ 



(15) 



where dnut denotes the partial derivative -rr-^, 

OXp 



4 = j/c e Z"' such that e{Q + k)nV 0| 



and where, for any p G M.'^, ij)p is the solution (unique up to the addition of a constant) to 



-div [AperVipp] = div [Bper {p + Vtf^)] 

ipp is Q -periodic, 



(16) 



and Xp is the unique solution to 



-div[AperVXp] = div [igBperip + Vw")] 

XpeLUR''), Vxpe{L\R''))\ 
lim Xp{x) = 0. 

V Ixl— >-oo 



(17) 



We assume that G W^'°°{V) andu{ G W^'°°{V). Then 



E 



(18) 



where C is a constant independent of e and rj. 



We wish to point out that the assumption Ug G W'^'°°{V) (and subsequently u\ G W'^'°°{V)) is a standard 
assumption when proving convergence rates of two-scale expansions (see e.g. [T, Theorem 2.1] and p. 28]). 
Note that, in view of ([2]), this assumption implies that the right hand side / in ([U]) belongs to L°°{T>). We also 
note that vf^ is not uniquely defined, since Wp and tpp are only defined up to an additive constant. However, 
adding a constant to any of these functions does not change the order of convergence in (jl8p with respect to e 
and r], but only the constant C. Choosing the best constants in Wp and ipp is hence irrelevant here, although 
it is an important matter from the practical viewpoint. Lastly, the existence and uniqueness of a function Xp 
satisfying ()17p is shown in Lemma [TOl below, in dimension d > 1. In dimension d = 1, the boundary conditions 
of (jl7p need to be modified for this problem to have a solution. The one-dimensional version of Theorem [2] is 
as follows: 
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Theorem 3. Assume that the dimension d is equal to one. Let he the solution to Q in the domain P with 
f G Lp'iV), and assume that satisfies (f9|)- (fT0|) . Let vf^ be defined by (fT5|) . where the definition pT|) of the 
function x is replaced by 



n' 



loc\ 



[l(0,l)-Sper 



;i + (u'0)')]' in 



where solves (H]). Th 



en 



E 



\u„ 



E 



< C (e + r/V^ + r?2) 

< C (e + r/V^ + 7?2) 



(19) 



(20) 



where C is a constant independent of e and r]. 



Note that, in dimension d = 1, we do not need to assume pT]) and (fT2]) . In dimensions d > 1, as pointed out 
above, these assumptions are used to have that the correctors Wp and tpp, solutions to dH) and (fT6]) respectively, 

In dimension d = 1, the coercivity assumption on Aper and the boundedness 



l,oo 



both belong to W 

assumption on -Bper are enough to show that and both belong to W^'°°(M.). Likewise, when d > 1, we 
assumed that Uq G W'^'°°{'D) and G 1^^'°°(2?) (which implies that / G L°°{'D)). When d= 1, the assumption 
/ G L'^(V) is enough. 

On another note, we notice that x is now only defined up to an additive constant. Again, changing x by a 
constant does not change the order of convergence in (fTn|) - (pU|) with respect to e and r], but only changes the 
constant C. 

In addition to its theoretical interest. Theorem [2] has also interesting numerical counterparts. Indeed, to 
compute vfj, one needs to solve problems set on a bounded domain (either with Dirichlet or periodic boundary 
conditions), and to solve for Xpi solution to the problem ()17p . set on the entire space. However, the right hand 
side in (fT7|) is the divergence of a compactly supported function, and we will see that Xpi^) quickly vanishes 
when X is sufficiently large (see Lemma [TOl below) . Hence, in practice, it is possible to approximate (|17p by 
using Dirichlet boundary conditions on a domain of limited a size. 

The proof of Theorem [2] consists of two steps. The first one is to expand with respect to r]. This 
is performed in Section [2] below (see Lemma H]). Each term of the expansion of is found to be the unique 
solution of a partial differential equation with a deterministic, highly oscillating coefficient, to which is associated 
a homogenized equation. The second step of the proof consists in successively estimating, for each of the terms 
of the expansion in rj, the rate of convergence of their two scale expansion in e. Corresponding results are stated 
in Section [3] (and proved in Section [5]). Collecting these results, we are then in position to prove our main result, 
Theorem [2] (see Section HI where we also prove Theorem [3|) . 



2 Expansion in powers of r] 

In this section, we expand the solution to ^ with respect to r/. 



Lemma 4. Let be the solution to 
follows: 



where Uq is solution to the deterministic problem ([T]), is solution to 



Under the assumption ([9]), it can be expanded in powers of r] as 

(21) 



-div 



Aper{-]Vul{-,u;) 



div 



Ai [-,u)Vul 



in P, 
on &D, 



(22) 
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and is solution to 



-div 



div 



in T>, 
on dD. 



(23) 



I r^(-,u;) = 
In addition, we have, almost surely, 

ll'^ol^M^') - ^' \\uii-^^)\\m{V) < C, < C, 

where C is a deterministic constant independent of e and rj. 

Proof. The relation (j21|) is a simple consequence of the linearity of the considered equation. The bounds 
follow from the uniform ellipticity of the matrices Arj and viper; and the boundedness of ^i. □ 

For the sequel, it is useful to further decompose uf in a deterministic part and a stochastic part of vanishing 
expectation. 



Lemma 5. Under assumptions ([9])- ()10p . the solution uf to ()22p writes 



where is the unique solution to 



-div 



Aper ( - ) Vu\ 



div 



Bper ( - ) VUq 



in D, 
on dD, 



and is the unique solution to 



-div 



div 







in T>, 
on dT). 



In addition, there exists C, independent of e, such that 



E 



< C. 



Assume furthermore that (fTTj) holds, and that 

f E L'^iV) for some q > d. 
Then there exists C, independent of k and e, such that 

W^IWl^CD) < Ce. 

Proof We note that, if A; G Z"* is such that e{Q + k)nl) = 9, then ^ writes 



-div 



A 



per 



V0I 



in V, 
on dV, 



(24) 



(25) 



(26) 



(27) 



(28) 



(29) 



the solution of which is obviously (pf. = 0. The sum in (|24p hence only contains a finite number of terms, and 
the proof of the decomposition ()24p goes by linearity of the equation (j22p . Note however that the number of 
terms in ()24p depends on e, and diverges when e — )• 0. 



We now prove the bound (|27p . Using that Ap^r is coercive, we infer from ()26p that 



< jJVcPlflQ+k [-) Bp,r [-) Vuf, 

< \\Bper\\L'^\\4'k\\H'^iV)\\Uo\\m(£(Q+k)), 

where a > is some constant that only depends on the coercivity constant of Ap^r and the Poincare constant 
of the domain D. Thus 



le ||2 



Using that (/)| = as soon as e{Q + k) CiV = 



'per||L°o ll'"0ll/i-i(e(Q+fc))- 



£||2 



we obtain 

-2| 



,e||2 



We deduce from that bound and the assumption that the random variables Xk are i.i.d. that 

2 



E 



Var(Xo) 5] ||</.|| 



ffi(X)) — Var(Xo)a ||-Bper 1 1 1°° 11^^011/^1(1)) — 



where C is independent of e (we have used (p3|) to bound Ug). We thus have shown ([27|) . 

We finally turn to the proof of ([29]) . Let us define (pl.{x) = (j)%{ex) on V/e. In view of ([26]) . we see that (pi. 
solves 



f -div 


Aper^ (pk 


= ediv [lQ+fcSperVMg(e-)] 


in P/e, 




D 




on d{V/e 



Introduce now the Green function Ti,{x,y) associated to the operator L = — div [Ape^V-] on the domain V/e, 
with homogeneous Dirichlet boundary conditions. We recall that r^(x,y) := r£(y,x) is the Green function as- 
sociated to the adjoint operator = —div [^pg^V-] on the domain T>/e, with homogeneous Dirichlet boundary 
conditions (a proof of this fact is given in [16t Theorem 1.3] and \\.2\ Theorem 1] in the case li > 3, and this 
proof carries over to the case d = 2). Consequently, we have Ti;{x,y) = as soon as x or y belongs to the 
boundary d(V/e) We can thus write 



e / ^e(x,y) divy[lQ+k{y)Bper{y)Vul(ey)] dy 

Iv/e 



-e I V yVe{x,y)Bper{y)Vul{ey) dy. 

'Q+k 



Hence, for any x £ V, we have 



^ X, 



e(Q+k) 



VyT, l^,^') Bper (|) Vu'o{y)dy. 



8 



Using the fact (see [3l Proposition 8]) that, under assumption ()lip . the Green function on the domain T) je 
satisfies 

G Vje, yy € V/e, |V,r,(x, y)| + |Vyre(x, y)\ < _ y^d-i (30) 
for a constant C independent of e, we have 

|(/.|(:c)| < C||5per||L-||Vug|Uoo(^) / ^ ^^^^2/- (31) 

Je(Q+fc) F ~ y\ 

We will show in the sequel that (|28p implies that there exists C such that, for all e, 

\WuI\\l^(V) < C. (32) 

We are thus left with bounding the integral in ()3ip . To this aim, we distinguish two cases. If |a; — ek\ < e, then 
there exists a constant than only depends on the dimension such that e{Q + k) C B{x,P(i£) (for instance, in 
dimension d = 2, p2 = 1 + \/2/2). We then have 



L 



I f 1 

-rdy < / rr~rdy < Ce, C independent of e. (33) 



(Q+fc) \x -yV ^ ' JB{x,pae) \^ - y 



Otherwise, if \x — ek\ > e, then any y E e{Q + k) satisfies p^e < \x — y\ for a constant p^ that only depends on 
d (in dimension d = 2, p2 = 1 — \/2/2). For those x, we have 

/ 1 rr~r — 7= — ^rr~r I dy < Ce, C independent of e. (34) 

JeiQ+k) F - y\ {Pd^r JeiQ+k) 

Thus, collecting ([31]), ([32]), ([33]) and ([Ml), we obtain that 

11*^1 IIl°°{d) — ^ independent of e. 

Proving ()29p therefore amounts to now proving (j32p . Again using the Green function r£(x,7/) associated to the 
operator L = — div [j4perV-] on the domain "D/e, with homogeneous Dirichlet boundary conditions, we write 

Vuf^ix) = e'-'' [ VxrJ-,y) f{y)dy. 

Using the bound (f30p . we deduce that there exists C independent of e such that 

VxeP, \Vuf,{x)\<C j^j^^^^dy. 

In view of assumption (j28p . we have / G L'^iV) for some q > d. Using Holder inequality, we write 

VxGP, |V^/^(x)| < C||/||i,(^) 



1 



Id-l 



1 1 

- + — = !■ 



(©) 



The function y i— )• |x — y]"^"'^ belongs to LP(T>) for any p < d/{d — 1). Since q > d, we have < d/{d — 1), and 
the norm in L"? of y i— >• |x — y|^~'^ is independent of x. The above estimate thus yields ()32p . This concludes the 
proof of Lemma [5l □ 
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3 Two scale expansions in powers of e 

Collecting (j2T]) and (pi]) , we have obtained that 



?x^(x,a;) = ul{x) + rj 



+ ri\'^{x,uj), (35) 



where is bounded in H^(T>) uniformly in e, rj and oj (see ()23p ). 

We now consider successively each term of the above series and show a rate of convergence on the difference 
between Uq, u\ and (/)| and their respective two-scale expansions. For clarity, the proofs of our results are 
postponed until Section [5l 

We start by Uq solution to ([1]). Note that this problem is a classical periodic homogenization problem, 
the limit of which, when e — )• 0, is well-known. The following result, giving a rate of convergence of Uq to its 
homogenized limit, is also classical (see e.g. [17, p. 28]). 

Proposition 6. Let Uq and Uq he the solution to ([1]) and respectively. For any p € M*^, we assume that the 
solution Wp to dH) satisfies G W^'°°{R^). We also assume that Uq G W'^'°°{'D). We then have 



d 

^0 = + ^ E < (- ) ^^"0 + eOl, (36) 



i=l 



where 6q satisfies 



for a constant C independent of e. 



\ee'o\\HHv) < C^e (37) 



We recall that, under assumption ([TT]) . we indeed have that uP^ E P^^'°°(R'^) for any p G M*^ (see e.g. [15} 
Theorem 8.22 and Corollary 8.36]). 

We now turn to lif solution to (j25p . This problem is not a classical homogenization problem, since its 
right-hand side also varies at the scale e, and only weakly converges in H^^{T)) when e — t- 0. We first proceed 
formally, using the two-scale ansatz approach, to identify the homogenized equation. We next state a precise 
homogenization result, and finally evaluate the rate of convergence of the two scale expansion. 

To derive formally the homogenized equation associated to ([25]) . we make the classical two-scale ansatz 



X \ -| / X \ 9 o / X 

i^{x) = Ui [x, — j + eul \ x, —j+eui yx, — 



where each term of the above expansion is assumed to be periodic with respect to the second variable. Inserting 
this ansatz in ()25p and using the two scale expansion ()36p of Uq (where we neglect the remainder eOq), we can 
easily derive a hierarchy of equations. We deduce from the equation of order e~'^ that is independent of its 
second variable: u\{x^y) = u\{x). The equation of order reads 



-divy [Aperiy) {V^u{{x) +Vyul{x,y))] = ^^diUoi^) divy [Bperiy) {ci + Vyu;° (y))] . 

1=1 

Using the functions Wp and t/^p defined by (H]) and (fT6]) . we thus see that 

d 

u\{x,y)=T{x) + Y,^iul{x)iPeM + ^^ul{x)wl{y), (38) 



i=l 
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where r is an undetermined function that only depends on x. We are now in position to use the equation of 
order e", which reads (recah we have neglected the remainder eO^ in (j36p ) 



diVa; [Aper{y) xu\{x) + Vyu\{x,y))] - divy [Aper{y) {yxu\{x,y) + Vyu\{x,y))] 



diVa; [Bper{y) [^i + VyWi" (y)) Sj-Ug] + divj^ 



i=l 



^wl^{y)BpeT{y)^ xdiul{x) 



i=l 



We close the hierarchy by integrating the above equation over the variable y £ Q, using that y i— t- nf(x,y) is 
Q-periodic. Using ([38]) and the expression ([3]), we then obtain that satisfies 



-div [ylpg^Vn|] = div 



BVu, 



= 



in V, 
on dV, 



with 



(39) 
(40) 



yi<i,j <d, B,j = / ejAperViJe, + / ejBperiej + Vw^). 
JQ JQ 

Mutiplying ([TU]) (for p = ej) by tOg. and integrating over Q, we find that 

VI <i,j<d, [ {VwlfAperVi^e, = " / {Vw^f Bp^ricj + Vw^). 
JQ JQ 

Inserting this relation in ()40p . we deduce that the matrix B is equal to the matrix B defined by ()13p . We hence 
deduce from ([39]) that indeed satisfies ([H]). 

These formal computations are formalized in a rigorous way in the following Propositions: 



Proposition 7. Assume that, for any p G M"^, the corrector Wp solution to @ satisfies Wp G W^'°°{R.'^), and 
that the solution Uq to ([2]) satisfies Uq G W^'°°{'D). Then the function uf solution to ()25p converges, weakly in 
H^{'D) and strongly in L^{'D), to the unique solution to (fT4 



The regularity assumptions on Wp and Uq ensure that Vuq in the right-hand side of ([25]) can be controlled 
in the appropriate norm. 

Proposition 8. Letu\ he the solution to ()25p . he the solution to ()14p and he the solution to ([2]). For any 

p G M'^, let he the solution to (jH) and ipp he the solution to (|16p . 
Introduce vf defined by 

d 

vl = u*^+eY^ ( - j diu{ + Ve, (-) diu^j , 



1=1 



and assume that Uq G W'^^'^{V), G W'^''^{V), and that, for any p G W^, we have G VF^'~(M'^) and 



1,CXD 



We then have 



\Wi-vI\\hHv) < C\/e 



for a constant C independent of e. 

Again, under assumptions pT]) and we have Wp G W^'°°{R'^) and i/jp G 
e.g. [la Theorem 8.22 and Corollary 8.36]). 

We finally turn to 0| solution to dJSDi namely 

—div 



for any p G M'^ (see 



Ap,r[-)V<Pl\=dW[cl] inV, 
on dV, 
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with 



- B 



per 



Assume momentarily that the sequence Vuq is bounded in L°°{'D) (we have proved such a bound above, see ([32]) . 
under the strong assumptions (|28p and (jlip l. Then, for any k G Z*^, c| converges to in L^(I?). Using the 
coercivity of Ap^r, this imphes that (p^. converges to in H^{T)). We thus have the following result, which will 
be rigourously proved in Section [S] below: 

Proposition 9. Let (/)| be the solution to (pUj) . and let Uq and Wp be the solutions to ([2]) and (H]). Assume that 
u* e W^'°°(V) and that, for any p G W^, we have w° G PF^'°°(M'^). Then cpl converges to in H^{V). 

To describe more precisely the behavior of (p^., we need to introduce the auxilliary function Xp defined 
by (mi below. Recah first that Q = (-1/2, 1/2)"^. Following the same arguments as in [21 Lemma 4], we have 
the following result, which will be useful in the sequel. 



Lemma 10. For any p G W^, the problem 



-div[AperVXp] = div [iqBperip + V 

XpeLUW''), VxpG(L2(M'^))^ 



(41) 



has a solution which is unique up to the addition of a constant. In addition, under assumption (jlip . there exists 
a solution of ()4ip and a constant C > such that 



Vx G 



(J 

with \x\ > 1, |VXp| < -j — rj, 



Vx 



C 



+ F 



(42) 
(43) 



In the sequel, we will always refer to that particular solution of (|4ip . 



We are now in position to make precise the behavior of i?i>| in the norm. Let us first argue formally. 
Introduce the matrix = Iq+kBper- Using the periodicity of Ap^r, Bp^r and w^, and after changing variables, 
we recast (HTll as 



-div 



A 



per 



- Vxp - - 



div 



In turn, the problem (|26p reads 



-div 



A. 



per 



div 



Ek (-) Vug « ^div Ek (-) diu^ (e, + V< ( - 



1=1 



where we have used the expansion ()36p of Uq (in which we have only kept the highest order terms). Assuming that, 

d 



in the above equation, x and x/e are independent variables, we thus see that V(/)|(x) diUQ{x) Vxe^ f k), 



i=l 



and thus 4>ki^) ^ ^^^'5jUo(x) Xet y kj. These formal manipulations motivate the following result, the rig- 

orous proof of which is postponed until Section (5) 

Proposition 11. Let (j)^^ be the solution to ()26p and Xei be the solution to (j4ip . for 1 <i < d. Introduce 

le = \k e 'L'^ such that e{Q + k) r\V ^ $^ , Card{I,) ^ e''^, (44) 
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and 



where Uq is solution to Assume that Uq £ W'^'°°(T>), and that (|lip holds. We then have 



where C is a constant independent of e. 



4 Proofs of Theorems [2] and [3] 

Proof of Theorem\^ We have shown above (see (|35p ) that 



u. 



^(x,a;) = ulix) + 7?E(Xo)nf (x) + r? ^(X,,(a;) - E(Xo))</.|(x) + r^\'^{x,u), 



where the set is defined by (j44p (recall that (/)| = whenever /c G Z'^ is such that k ^ I^). Using the fact that 
Xk are i.i.d. scalar random variables, we have 



E 



< C [Dl + Dl + Dl + Z)i] , 



where 



p=i 



Dl = 7?|E(Xo) 



^! - - (< (-) + V'e, (-) dpUl 

p=l 



I^s = 77v'Var(Xo: 



E 



r?%/E 



'»7llHi(D) 



We have shown in Propositions [6l [8l and [TT] that ^ C^/e, Di < Cr]\f£ and D2 < Cr\^Je ln(l/e) respectively, 
for a constant C independent of e and 77 (note that all assumptions of these propositions are satisfied since, in 
view of (HI]) and ([l2|), we have w° € VF^'°°(M'^) and Vp G VF^'°°(M'^) for any p G M"^). Next, using Lemmagl we 
see that -D3 < Cry^ for a constant C independent of e and ry. This concludes the proof of □ 



Proof of Theorem\^ To fix the idea, we choose T> = (0, 1). We again argue on the basis of (|35p . Tedious but 
straightforward computations show that, in dimension one, the estimates of Propositions [6l [HI and [TT] read 



dx 



< Ce, 



L2(0,l) 



dul 


dUf 


dx 


dx 



< Ce, 



L2(o,l) 





_ M 


dx 


dx 



< Ce. 



L2(0,l) 
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We thus obtain 







duf^ 




2 


E 












dx 


dx 





< C7 (e + + r/2) . 



We next write that, almost surely, 



L°°(0,1) 



< 



du^ dv^ ^ 



L2(0,l) 



+ |n^(0,w) -^;^(0,L^)| 



Using that u^(0,a;) = Uq{Q) = il^(O) = and that and ip belong to L°°(M), we obtain that 



4(0,a;) -t;^(0,a;)| < Ce + Cer] 



S)'(0)^(X,(a;)-E(Xo))x(-fc) 



hence, using that x ^ 
E 



we have 



(45) 



(46) 



Collecting this result with ()45|1 and (|46|1 yields the bound (pOj) . Likewise, collecting (|20|) and (|45]) . we obtain 
the bound p^ . This concludes the proof of Theorem [3l □ 



5 Proofs of the two scale expansions 

We collect in this section the proofs of the results stated in Section [3l The following technical result, already 
present in [n\ p. 27], and that we recall here for the sake of completeness, will be useful. 



Lemma 12. Let V be a bounded open set ofR'^. Consider Z G (L^ 



loc 



a Q -periodic vector field such that 



div (Z) = and Z = 0. 

Jq 



Then, for any v G W '°°{T>), we have 



div 














H 



< CeWVvl 



where C is a constant independent of e and v. 
Note that, as Z is divergence free, we have 

div 



Z[-]v 



Z(-\-Vv. 



Since Z is Q-periodic, this quantity converges weakly in L^iV) to {Z) -Vv = 0, as the average of Z vanishes. 
The above result hence shows that, in the H~^{'D) norm, the above quantity vanishes at the rate e. 

Proof. In view of the assumptions of Z, there exists (see [171 p. 6]) a skew symmetric matrix J such that, 



VI < i < d, Zj = 



i=l 



dxi 
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and 



yi<i,j<d, Jij e hI^cO^'^): Jjj is Q-periodic, / Jij = 0. 

Jq 



The j-ih. coordinate of the vector Z [-] v reads 



Z(-)i,(i) 



1=1 



d 



V[X) — £ 



i=l 



eBj{x) — eBj{x), 



where 



S,(x) = ^J.,(^)|^(x) and B,{x) = j2^^{J^,{^)H^))■ 



.£/ dxi 

1=1 i=l 



The vector B{x) is divergence free as J is skew symmetric. For any (p G Hq[T)), we thus have 



div 



-e (div [B] , , 

e / B-V(t> 
Jv 

d 



^ 1.7 — 1 



hence 



div 



Z[-\v 



< e||Vv||Loop)||(/)||^i(P) ^ (- 

d 

< e||Vv||ioo(X5)||(?!)||//i(x)) ^ ll<^ijllL2(Q) • 

As the above bound holds for any (/> € ifQ(P), we deduce that there exists C such that, for any v S W^'°°{T>) 
and any e, we have 



div 














H 



This concludes the proof. 



□ 



5.1 Two scale expansion of u\ 

In this section, we prove Propositions [7] and [8l 

Proof of Proposition This homogenization result is proved using the method of oscillating test functions 
[25 1 . The variational formulation of (1251) reads 



(47) 
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where, for any u and v in Hq{T>), 

Ae{u, V) = (Vvf Aper (-) Vu and Ls{v) = {Vvf Bper (-) Vug. 

Using the coercivity of Aper, the boundedness of Bper and ([23|) . and taking v = uf as a function test in (f^7|l . we 
obtain that uf is bounded in H^(p). Thus, using the ReUich Theorem, we deduce that there exists G Hq(T>) 
such that, up to the extraction of a subsequence, 

ul converges to u^, weakly in Hq{T>) and strongly in L/^iJ)). 

For any function (/? G C^(2?), define the test function 



i=l 



which obviously belongs to Hq{T>). In view of (j47p . we have 

A(«!,^^^) = -^.(^^'). (48) 
We now expand both sides of ([^8]) in powers of e: 

Ae{u\,v') = Al{u\,if) + eAl{ul^), (49) 
L,{v^) = Ll{^) + eL\{^), (50) 



where 



^,i(n5, ^) = f (;) (^^^^)^ ( - ) 

i=l 



„ CI 



We now successively study the limit of these four quantities as e — )• 0. Using ([23]) . the fact that tfg. E 
that lif is bounded in H'^{V) and the boundedness of Ap^r and -Bpen we obtain 



l,oo (TBd\ 



\A'^{u{,'f)\<C and |L^((^)| < C, C independent of e. (51) 
We now turn to L^. Using the two scale expansion (j36p of Ug, we see that 

LO(^) = (v.) + LOi(^) + Lf{^), (52) 
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where 

Lf{^) = E / (e. + V< {-)YB,,r (-) (e, + Vw^ (-)) dX„ 
Lf{^) = ^Y. j (e^ + V<(-))^i?per (-)<(-) (V9,n*) 

Using dSZD, w° G W^'°°{R'^) and G VF2'°°(P), we obtain that 

\Lf{ip)\<C\\e9'o\\HHv)<CV^ and \Lf{ip)\<Ce, (53) 
where C is a constant independent of e. Turning to L^*^, we see, using that Bp^r and Wg. are Q-periodic, that 

(e, + Vu;0 (^-) ) ^ (-) [ej + Vu;°^, (-) ) ^ weakly-* in L°°, 
where -B is defined by ()13p . Thus 

^ / (V(^)^WnS ase^O. (54) 
Jv 

Collecting ([52]), ([53]) and we obtain that 

L°((^) ^ / {VipfBVu^ ase-^0. (55) 



We next turn to A^. Using that div viper j (^i + VtOg^ (~ ^ ^ ~ ^ symmetric, we obtain 



that 



d 



(56) 



Recall now that uf — t- strongly in L {T>) and that, as Ap^r and z/;^. are Q-periodic, we have 
Aper (-) (ci + Vu;° (~) ) ^ ^per (cj + Vw° ) = vlp^^ei weakly-* in L°°, 
where j4pg^ is defined by ([3]). We thus deduce from ([56]) that 

d „ 

Al{ul,^p) ^ I u^{Vdnp)'^A* a ase->0. 
i=i -^^ 

Collecting ([18|) . (09]), (|50|) . ([5T]) . the above limit and ([55]) . we obtain that satisfies 

-E / ^t(Va,^)'^4g,e, = - / {V^fBVu*, 
f^^Jv Jv 

for any ip £ Cq^{'D). This shows that solves (which has a unique solution) and thus concludes the proof 
of Proposition [71 □ 
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Proof of Proposition\^ The proof mostly goes by using the coercivity of Ap^r and showing that, in some ap- 
propriate norm, — div [Aper(Vuf — Vuf )] is small. However, a technical difficulty comes from the fact that 
Vi ^ HqIV), as it does not vanish on dV. A preliminary step (Step 1 below) thus consists in approximating 
by a function (namely gf defined by (j57p below) that is equal to vf away from the boundary dT), but vanishes 
on the boundary. Step 2 consists in estimating the difference uf — gf- 

Step 1: Truncation oi vf 

Let us define G C^{T>) such that < Ti:{x) < 1 for all x £ D, Tei^) = 1 when dist(9P,x) > e and 
^ll^''"e||L°°(i') ^ C*, where C is a constant independent of e. We denote by C P the set of M*^ defined by 

:= {x G P such that dist(5P,2;) > e} 

and we note that 

\V\V,\ < Ce. 

Introduce now gi G Hq{T>) defined by 

d 

gl = u^^ + eT,Y, (< (-) + V'e, (-) diui) , (57) 



where Ui is the solution to ([H]), Ug is the solution to ([2]), and tt^g. and ijjei are solutions (with p = ei) to @ 
and (|16p . respectively. Note that g\ = Uf except in a neighboorhood of 92?. In the sequel, we estimate v\ — gf- 
In the next Step, we estimate gf — uf. 
By definition, 

Vvl-Vgl = el-el + eel (58) 



where 



a 

el = (1 - r,)j;(v<(-)d,n^ + We, (-) 5,^x5) 
1=1 

d 

1=1 

d 

e| = (l-r,)^(<(-)v(a,nt) + ^e,(-)v(a,nS)). 



1=1 

i=l 

We now bound from above successively the norm of el, ef and Cq. First, as Uq G H^^'°°(2?), G VF^'°°(P), 
V'e, G T^^'°°(M'^), u;0 G T^i'°°(M'^) and < < 1, we have 

||e||||2(x)) ^ C*, C independent of e. (59) 

The same arguments lead to 

2 



|6iIIl2(d) 



.1=1 

■el 



< Ce, (60) 
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for a constant C independent of e. We next write 

d 



i=l 



< Ce. 



L°°{V) 



Collecting dSHD, ([59]), ([601) and ([SI]), we have 



\Vv\ — Vgf ||^2(-p) < Ce, C independent of e. 



Observing that 



(61) 



a 



i=l 



we obtain that 



1 ~ 9i\\h^(v) ^ C-v/e, C independent of e. 



(62) 



Step 2: We next turn to estimating -uf — Qi- Using that Ap^r is coercive and the fact that u\— ^ Hq{'D), we 
have 



o^\\Ui — 9i\\h^(V) 



< 



[ (VI 
Jv 



v\ - VglY Aper - (Vnf - V^f ) 



< (Vuf - VglY Aper (-) (Vnf - Vv\) + (Vnf - Vg{)^ Ap^r (-) (Vl^f - V^f ) 



< \Wi-9i\\m(v) 



div 



Aperi-] (VHf - VTJf) 



■per\\L°°\\vl — gi\\m{V) )(^3) 



where the constant a > only depends on the coercivity constant of ^per and the Poincare constant of the 



domain D. In the sequel, we bound from above 
By definition of , we have 



div 



per 

Aper (-) {Vul - Vvf] 



v\ = v{ + v{, 



with 



\ + e^wl(-\diu{ and vl = e^il^e, [-j diU^. 



i=l 



1=1 



Using the equation ([25]) on u\ and the relation ([H]) between u\ and Uq, we compute 



div 



Aperi-] (Vuf - Vuf) 



div 



A, 



-per 



Vv{ - A%,Vu\ 



+ div 



A^ 



per 



Vvl + Bper ( - ) Vug - BVUq 



Dq + Di+ eD2 



(64) 



where 



Do 



j=i 



4 



ei + Vu;°. 



4* p- 



1 ' 



A, 



per \ - J VV'ei 



per 



Si + Vw" 



I?2 



div 



i?pe.(-) V^gJ + J^d 



IV 



per 



Vdiul + V'e, (-) V9i<) + Sper (-) < (-) V^iU* 
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We now bound from above these three quantities. As Aper and -Bper are bounded, we see that 

d 

\\D2\\h-^V) < C\\9o\\h1{V) + C'X] [ll^eJI-f'°°ll'"lll^^2(C) + \\^eML^\\Uo\\H^(^v) + 1^0 II ^2(7?)] : 

i=l 

from which we infer, in view of (j37p . that 

£111)2 ||_H--i(D) < C-v/e, C independent of e. (65) 

Let us now turn to Dq. Consider, for any 1 < i < d, the vector- valued function 

Z{y) = Aperiy) (e, + Vwliy)) - A^^^ei. 

We observe that Z e {Lf^^{R'^))'^ is diver gence free, Q-periodic and of vanishing mean. Since diu{ S W^'^{'D), 
we can use Lemma [T2l and we obtain 

||Do||i?-i(X') ^ Ce, C independent of e. (66) 
Turning now to Di, we hkewise consider, for any 1 < i < d, the vector- valued function 

Z{y) = Aper(2/)VV'e,(y) + 5per(y) (ci + V<(y)) - BCi. 

By construction, Z G {Lf^^(R'^))'^ is Q-periodic and divergence free, in view of the definition ()16p of ipet- In 
addition, the mean of Z vanishes. Indeed, for any 1 < j < d, using ()13p . ()16p . the symmetry of ^per and Q, 
we have 

Z-ej = / ejAperVV'e, + / cJSper- (e^ + Vii;° ) - / {ej + Vw^ fBper{ei + Vwl) 
JQ Jq Jq Jq 

= [ ejAper-VV'e, - / (Vu-I?. )^5pe, (e^ + Vu;° ) 
Jq Jq 

= / ejAperVtpe,+ / (Vu-^ . )^^per- V^'e, 
(VVeJ%er(ej+VO 



Q 

= 0. 

Since diU^ S Ty^'°°(P), we have that Z and SjUq satisfy the assumptions of Lemma [121 hence 

ll-^i||_f/-i(r') ^ C'e, C independent of e. (67) 
Collecting ([Ml), (ESI), dSSI) and dSZD, we have 



div 



where C is a constant independent of e. We now infer from and pSj) that 



an 



■E „£||2 



5lll_H-i(») - ^11^1 ~ 5'lll/i'i(D) V^) 



'1 - yiiiHi(») 
hence 

11^1 ~ 9i\\h'^{V) ^ Cy^) C independent of e. 

Step 3: Conclusion 

Collecting the above bound with ([52]) . we deduce that 

11^1 ~^i\\h^{v) ^ Cy^) C independent of e. 
We thus have proved the claimed bound, and this concludes the proof of Proposition [HI □ 
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5.2 Two-scale expansion of 0| 

In this section, we prove Propositions [91 and 1111 
Proof of Proposition O Introducing 

cUx) = iQ+fc (^) Bper (^) Vn^(x), 

the problem ([26|) writes 

/ -div \Aper (^) V^l] = div [c|] in P, 



on 9P. 



Multiplying this equation by cpf,, integrating over V, and using the coercivity of Ap^r, we obtain that there exists 
C independent of k and e such that 

m\\miV)<C\\cl\\L2(^T,y (69) 
Let us now show that c| converges to in Lp'iT)). Using the expansion ()36p . we write 

with 

d d 

= E^^^S (e^ + V< (-)) +e j;V(a.nSX (-) . 

i=l i=l 

Using the fact that wl G VF^'°°(]R'^) and G we see that is bounded in L°°{V). We next write 

ll^£||2 ^ II R ||2 f lY7o,e|2 

d , I l^v7fl£|2 



^e{Q+fc) 

< Ce"^ + C||e6'o||^l(X))• 
Using the bound (f37|l . we deduce that c| converges to in L^(P). In view of ([69]) . this implies that 0| converges 
to in Hq{'D). This concludes the proof. □ 



Proof of Proposition\ll[ As in the proof of Proposition [HI the proof falls in two steps. We first truncate v^. in a 
function vf. (defined by ()70p below) that vanishes on dD. We next estimate the difference between and (/>|. 

Step 1: Truncation of I;| 

Let us define G C^{T>) such that < Ts{x) < 1 for all x £ D, Tei-i^) = 1 when dist{dT>,x) > e and 
£||V''"£||l°°(X') ^ C*, where C is a constant independent of e. We introduce 

Ve := {x G P such that dist(aP, x) > e} , \V\ V^l ~ e, 
Je := |A:GZ'^suchthate(Q + A:)nP\Pe/0}, Card(Je) ~ e^"'^, 

and the function vf^ G Hq (D) defined by 

vl = eT,Y,Xe,{--k)d,u*, (70) 
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where Uq is solution to ([2]) and is solution to ()4ip . Note that vf, = except in the neighboorhood 
boundar 

Yl \\<Pk ~^k\\m{v)^ where, we recah (see (|MD), 



of the boundary of V. In the sequel, we estimate ^ — "wlH^i/x))' smd, in the next Step, we estimate 



4 = j/c G Z'^ such that e(Q + A;) n P / 0} , Card(4) ~ e" 



Recall also that, whenever k ^ Is, we have 0| = 0. 
By definition. 



V^)|-Viy| = eS'^-e^^ + e^•^ (71) 



where 



k,£ 



k,£ 



k,£ 



1=1 

e(l-Te)^Xe. (--k)vidiU*o). 



.=1 



We now bound from above successively the norm of 62'^, e^"^ and Cq'^. To this aim, the following computation 
will be useful: for any 1 < i < d, we have 



There exists p such that 

Ve, Vj e Je, Vfe G /e, Q + j-kc B{0, p/e). 



/(.ei 



We thus obtain that 



E/ 4(7-^)^E^'/ Xl<e[ xl, (72) 



^^j^JV\V, ^ ^.gj^ ■JB{0,p/e) JB{0,p/6) 

We next infer from ()43p that 



Jb(0,p e) Jb(0,p e) U + Wl ) Jl ^ 



'B{0,p/£) ' JB{0,p/e) (1 + \y 

where C is a constant independent of e and 



(73) 



Collecting (j72|) and ([73]) . we deduce that 



keh 
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We now bound eY- As G W'^'°°{V), and using ()75p . we have 



Ell fc,£||2 
F2 IIl2(D) 



j=i 



< Ce^||V\5lli^EE/^ 



fcG-fe i=l 



k e 

We next turn to e^^' . The same arguments and the fact that e||Vre||j;,=o < C lead to 



2 



,i=l 



i=i -^^^^^ 

< CeRde^ 



where we have again used (|75p . Turning to , we have, using Vxe^ £ (L^(]R'^)) 



fcG/s i=l keh 



j 









^ ^EE-'^E/ ivxej 

< C^EE^'ll^^eJli. 

< Ce. 



Collecting (HI]), ([76]), dZZD and ([78]), we deduce that 

E W'^^k - Vi^||li2(i,) <C{e + eRd,e + e'Rd,e) 



keie 



where C is a constant independent of e. Observing that 



E 

keh 



I'^fc ~ ^fcllL2(-p) — C'e^||Vt(.Q||^oo(x)) 



7 ^ 7 > / , '^e, 



we obtain that 



[1 + ln(l/e)] if d = 2 
Ceifd> 2, 



where C is a constant independent of e. 
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Step 2: We next turn to estimating ^ ~ ^Illi/ifX))' Using that Aper is coercive and the fact that — f| G 

fee/, 

Hq{V), we have 



Vvlf Ap,r (-) (V(/)| - Vvl) = dI^' + 



(80) 



where the constant a > only depends on the coercivity constant of Ap^r and the Poincare constant of the 
domain 2?, and where 



D 





k,£ 
1 



(V0I - Vi^D^^per (-) (VcAl - Vlffc) , 
- VI^I)^ ^per (-) - Vvl) . 



We successively bound Dq'^ and dJ^'^ from above. We begin with Dq'^. Observe that, in view of ([26]) . 



div 



Aper(-) {V<Pl-Vvl) 



div 



-KB. 



hfc l^^y -per 



+ ^div Zk[-)dpul 



p=i 

+ ediv 



where the vector-valued function is defined by 

Zk{y) = ^Q+k{y)Bper{y) Up + Vu'°^(y)j + Aper{v)V Xe^iV " k). 



Note that, in view of ()4ip . is a divergence free vector, hence div 



Zk[-] dpUQ 



thus rewrite D^'^ as 



where 



-^fc (^-j • VSpUQ. We can 



(81) 



00 



(V0I - V?J|)^ iQ+fc (-) Bper (-) Viig - ^ (ep + Vwl (-)) apu^ 



k,£ 
01 



02 



-ej] / {V4>l-VvlfAper(-)v{dp 



^o) Xep ( - - A:) . 



We successively bound these three quantities. Since Xcp ^ L°°(M'^) (see (fl3|) ) and Uq G VF"^'°°(P), we have 
||-u|||ioo(x)) < Ce. We also have that ||</>|||loo(2?) < Ce, in view of (j29|) (recah indeed that € VF^'°°(P) implies 
that / € L°°{'D), in view of ([2]); assumptions of Lemma [5] are thus satisfied). Using that Uq E W'^'°°{V), we 
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now bound from above D, 



k,e 
01 • 



l^ofl ^ '^\\4>1-vI\\l°-{v)\N'^u*o\\l°-(v) 



p=l 

d 

< CeY. 

p=i 

d 

< CsY. 

p=l 



LHV) 



\B, 



per||L°° / iQ+fc 



+ ||^per||L°° 



< 



JV/e-k 

|VXeJ+/ 
\JB(0.1) JB 



IVXeJ 



'-8(0,1) JB(0,p/e)\B(0,l) 

Using that Vxe^ G (L^(M'^))'^ (see LemmafTUD and the bound (|^ . we deduce that 

^l/e 



p=l 

< Ce'^+i [1 + ln(l/e)] , 



1+ ( l|VXeJL2(Kd)+C7^ 



where C is a constant independent of e. We thus get 



Y,\Dlt\<Ce[l + H^/e)]. 



We now turn to D^f ■ Using ()73p . we observe that, for any k ^ I;, 



Xe,[--k 



Xe 



I |Xe,r<eM \Xe,\' <Ce''Rd,e- 

Iv/e-k JB{Q,p/e) 



We thus can bound from above -Dg^' using that Uq £ W'^'-^{'D): 

d 

fee/e p=i fce/s 



< Ce /j;||V</.|-ViJ|| 



\ p=l 



2 

L2{75) 



where C is a constant independent of e. We next turn to D^^. Using the bound (|37|) on the two scale 
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of Uq, we have 



Y.\D^00\ < ll^perlU- IIV0I - V?y|||i2p) 



p=l 



L2(e(Q+fe)) 



< C /5^||V0|-V^ 



r£ii2 



a 



p=i 



L2(e(Q+A;)) 



< C /5^||V0|-V^ 



fcllL2(x)) 



< CV~ejY.\\Vc^l-Vvl\ 



Collecting (IMI), dSl, (IS31) and ([51]), we obtain that 



V kale 



kei. 



We now turn to D-j^'^. Using ()79p . we have 



^e||2 



fcG/g V fce/. V fee/, 



^livij|-viy|||2 



L2(^) 



keh 



Collecting ([HOD, ([85]) and dMD, we obtain 



(84) 



(85) 



(86) 



a 



rre||2 



'kWV^V) 



with, in view of ([7H) . i?^^^ = 1 + ln(l/e) if d = 2, and Rd^^ = 1 if d > 2. This implies that 

Y2 W^k ~ '^kW'n'^ {V) — Celn{l/e), C independent of e. 

k&h 

Collecting this bound with (j79p . we obtain the claimed bound. This concludes the proof of Proposition [TTJ □ 
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